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ABSTRACT 


The  irradiation  by  infrared  laser  light  of  an  atom  adsorbed  on  the  surface 
of  a  harmonic  crystal  is  considered.  The  dynamic  coupling  between  the  atom  (its 
motion)  and  the  substrate  degrees  of  freedom  (phonon  field)  in  the  presence  of  a 
confining  potential  well  (van  der  Waals  bond  between  atom  and  crystal)  gives 
rise  to  thermal  relaxation  of  this  adbond  configuration.  Both  the  atom  and  the 
substrate  are  assumed  to  be  transparent,  but  the  bond  is  allowed  to  have  non¬ 
vanishing  dipole-moment  matrix  elements,  which  couple  the  external  field  to  the 
adsorbate.  The  equation  of  motion  for  the  reduced  adbond  density  operator  is 
obtained  with  reservoir  theory,  and  the  relaxation  constants  are  expressed  in 
properties  of  the  crystal.  With  a  similar  method,  the  spectral  profile  for 
absorption  of  weak  radiation  is  derived.  Subsequently,  the  illumination  by  a 
strong  f inite-linewidth  laser  field  which  is  in  close  resonance  with  a  single 
transition  of  the  adbond  is  examined.  The  optical  Bloch  equations  in  operator 
form  are  derived  and  applied  to  study  the  process  of  laser  heating  of  the 
crystal.  It  is  pointed  out  how  this  mechanism  can  be  understood  as  resulting 
from  (multi )photon-phonon  conversion  reactions  which  are  mediated  by  the  adbond. 


I.  INTRODUCTION 


Applications  of  tunable  laser  sources  in  contemporary  experiments  in 
physics  and  chemistry  ca.  roughly  be  subdivided  into  two  categories.  In  the 
first  class  the  laser  is  used  as  a  diagnostic  tool  to  investigate  a  dynamical 
system  with  spectroscopic  methods.  Scanning  the  frequency  of  a  low-power  field 
which  is  incident  on  an  atom  or  molecule  reveals  the  resonances  of  the  system  as 
lines  in  the  absorption  profile.  Interaction  of  the  system  with  the  environment 
(collisions,  spontaneous  emission,  presence  of  boundaries,  etc.)  amounts  to 
relaxation  of  the  molecule,  which  in  turn  broadens  a  spectral  line  in  such  a  way 
that  the  width  of  the  line  is  proportional  to  the  inverse  relaxation  rate  of  the 
specific  transition.  More  subtle  properties  of  the  line  shape,  such  as  the  far- 
wing  decay,  carry  information  on  the  details  of  the  interaction,  since  the 
frequency  dependence  of  the  absorption  is  determined  by  the  Four ier-Laplace 
transform  of  the  time-evolution  operator  of  the  density  operator.*  Therefore, 
it  can  be  expected  that  a  profound  comprehension  of  dynamical  properties  can  be 
achieved  from  an  accurate  observation  of  spectral  profiles.  However,  this 
method  requires  an  elaborate  theory  which  disentangles  the  contributions  of  the 
various  mechanisms  to  the  line  shape. 

In  the  second  category  of  experiments,  a  radiation  field  is  applied  in 
order  to  modify  or  even  Induce  a  process.  An  incident  photon  can  supply  the 
necessary  activation  energy  for  a  reaction  which  is  very  unlikely  to  take  place 
without  a  field.  The  advantage  of  a  narrow-band  laser  is  that  by  tuning  the 
laser  frequency,  the  process  can  be  optimized,  in  contrast  to,  for  instance, 
thermal  excitation  of  a  species  (collisions).  Prime  examples  are  laser-induced 
dissociation  and  ionization  and  atomic  resonance  fluorescence.  In  more  rigorous 
situations,  a  strong  (pulsed)  laser  is  merely  used  as  a  heat  gun  (melting. 


vaporization),  and  the  color  and  polarization  of  the  field  are  not  of  crucial 
importance. ^ 

In  the  present  paper  we  consider  optical  features  of  atoms  which  are 
adsorbed  on  a  crystal.  An  atom  is  bounded  to  the  substrate  by  electromagnetic 
interaction  with  the  crystal  atoms  (van  der  Waals  energy),  which  is  effectively 
described  by  a  potential  well.  The  few  (-25)  vibrational  adbond  states  have 
transition  frequencies  (level  separations)  of  the  order  of  50  cm  *  up  to  500 
cm  *,  which  is  in  the  infrared  (IR)  region  of  the  optical  spectrum.  Both  the 
atom  and  the  crystal  are  supposed  to  be  transparent  for  IR  light,  as  is  for 
instance  the  case  for  noble-gas  atoms  on  potassium  chloride.  The  vibrational 
states,  however,  exhibit  an  optical  activity  in  the  sense  that  there  are  non¬ 
vanishing  dipole  matrix  elements  between  the  various  bound  states.  Irradiation 
with  an  IR  laser  will  then  amount  to  photon  absorption  by  the  bond,  and  hence 
the  spectral  profile  will  reflect  the  properties  of  the  binding  potential  and 
the  dynamical  interaction  with  the  substrate.  In  this  fashion  we  can  study 
transparent  crystals  with  spectroscopic  methods,  where  the  medium  for  transport 
of  the  information  is  furnished  by  the  adsorbed  atoms.  This  situation  is 
reminiscent  of  the  more  familiar  problem  of  collisional  redistribution  by 
neutrals  in  gas-phase  experiments.  Transparent  atoms  are  immersed  in  a  buffer 
gas,  and  photons  from  an  incident  field  are  absorbed  during  collisions  only. 
Accurate  measurements  of  line  shapes  then  allow  the  determination  of  interaction 
potentials.  The  only  distinction  is  that  collisional  redistribution  deals  with 
scattering  states,  whereas  the  spectroscopy  of  adsorbates  involves  bound  states. 
It  can  be  anticipated  that  the  established  methods  for  the  gas  phase  can  be 

converted  into  suitable  surface  equivalents,  although  experiments  along  this 

3,4 


line  are  still  rare. 


Efforts  in  the  second  category  of  the  application  of  lasers  in  this  systen 

are  far  more  common,  and  most  notably  studied  is  the  process  of  resonant 

desorption.^  ^  An  intense  IR  laser  is  tuned  into  resonance  with  the  transition 

frequency  between  a  low-lying  and  a  high- lying  vibrational  state.  Absorption  of 

a  photon  is  accompanied  by  a  transition  of  the  adbond  to  the  high-lying  state, 

which  is  close  to  the  continuum.  Then  the  crystal,  regarded  as  a  thermal  bath, 

has  only  to  provide  a  small  amount  of  energy  in  order  to  accomplish  the 

desorption  of  the  atom.  Without  the  driving  laser  this  process  would  not  occur, 

since  the  atomic  bond  has  to  be  excited  to  the  high-lying  state  by  thermal 

coupling  to  the  substrate.  In  this  way,  one  has  expected  to  be  able  to  make 

clean  surfaces,  without  heating  (as  in  thermal  desorption)  or  damaging  the 
8-13 

materials.  Another  idea  has  been  that  advantage  could  be  taken  of  the 

resonant  nature  of  the  process.  Different  adsorbed  species  would  be  desorbed 

selectively,  depending  on  the  laser  frequency,  which  would  provide  a  practical 

12-13 

method  for  isotope  separation,  or  less  ambitious,  separation  of  molecules. 
Although  the  sketched  process  will  obviously  happen  if  the  radiative  coupling  is 
sufficiently  strong,  it  does  not  necessarily  mean  that  the  efficiency  is  very 
high.  Recent  experiments  show  a  quantum  yield  of  about  1 7.  (ratio  of  desorbed 
atoms  to  absorbed  photons),  which  conversely  implies  that  almost  all  radiation 
is  converted  into  thermal  energy  of  the  substrate.  Consequently,  the  solid 
heats  up  very  fast,  and  this  is  precisely  what  one  tries  to  avoid.  Since  the 
thermal  coupling  is  inevitable,  the  conclusion  is  that  resonant  desorption  is 
not  a  very  promising  technique.  Furthermore,  the  selectivity  has  turned  out  to 
be  very  poor,  if  present  at  all,  which  is  probably  due  to  a  rapid  energy 

14 

exchange  between  different  adsorbates  preceeding  the  desorption.  From  a 
different  point  of  view,  however,  this  mechanism  could  be  applied  for  laser 
heating  of  a  transparent  crystal.  By  manipulating  the  spatial  variations  of  the 


5 


laser  intensity,  this  process  can  be  used  to  maintain  well-controlled 
temperature  gradients  along  a  gas-solid  interface. 

II.  KINETIC  PHONON  COUPLING 

An  atom  with  mass  m  is  physisorbed  on  the  surface  (xy-plane)  of  a  crystal, 
and  its  motion  is  assumed  to  be  mainly  perpendicular  to  the  surface  (z- 
direction).  The  interaction  potential  with  the  solid  is  subdivided  into  two 
parts,  and  V^,  where  represents  the  potential  well  resulting  from  the 
coupling  to  the  closest  surface  atom,  and  is  an  effective  repulsive 
potential,  which  is  brought  about  by  the  remainder  of  the  crystal  and  confines 
the  atom  to  the  region  z  >  0.  The  thermal-equilibrium  position  of  the  nearest 
surface  atom,  with  mass  M,  is  taken  as  the  origin  of  the  coordinate  system. 

Then  the  positions  of  m  and  M  are  denoted  by  zez  and  u,  respectively,  and 
depends  on  | ze^  -  u|  only,  whereas  depends  on  z.  Subsequently,  it  can  be 

O  O 

asserted  that  u  (-  0.1  A)  is  much  smaller  than  z  (>  1  A),  which  allows  a  Taylor 
expansion  of  around  u  =  0.  Hence  we  can  write 


dV 

Vl(|zez  -  u|)  -  VL(z)  -  u-ez  +  . . 


(2.1) 


and  then  omit  the  ellipsis.  In  physical  terms  this  means  that  we  neglect 
multiphonon  processes  in  comparison  with  single-phonon  transitions,  as  will 
become  clear  in  due  course. 

Next  we  notice  that  the  expansion  (2.1)  separates  the  coupling  with  the 
lattice  vibrations  (last  term)  from  the  interaction  V^(z)  with  the  static 
crystal,  which  equals  the  interaction  for  u  =  0.  A  convenient  basis  set  for  the 
adbond  wave  function  is  therefore  provided  by  the  eigenstates  of  the  adbond 


Hamiltonian  H  ,  defined  as 
a 


where  the  first  term  is  the  kinetic  energy.  Once  a  potential  is  prescribed,  the 
eigenvalues  hw.  and  the  eigenfunctions  | i>  of  the  H  can  be  evaluated  with 

standard  methods.  A  common  choice  for  the  combination  V^(z)  +  V ^ (z)  is  a  Morse 

15-19  20 

potential  or  a  truncated  harmonic  potential.  We  shall  not  refer  to  a 

specific  potential  in  this  paper,  but  only  use  the  eigenvalue  equation 

H  |i>  -  hui.  |i>  (2.3) 

and  assume  the  eigenstates  to  be  non-degenerate  (as  for  a  Morse  potential). 

Then  in  turn  we  can  represent  Ha  with  respect  to  its  own  eigenstates  as 

H  =  V  rtw.P.  (2.4) 

a  L  1  i 

i 

with  =  |i><i|  the  projector  on  the  i-th  adbond  state. 

As  a  model  for  the  crystal  we  adopt  a  harmonic-lattice  representation,  for 

21 

which  the  Hamiltonian  reads 

Hp  -  J  lMk)aJsa  (2.5) 

ks 


in  terms  of  the  annihilation  and  creation  operators  for  phonons  (bosons)  in  the 
node  ks.  The  summation  runs  over  the  wave  vectors  k  and  branches  s  which  are 
supported  by  the  crystal,  and  the  dispersion  relation  w  *  w(k)  is  taken  to  be 
independent  of  s  and  of  the  direction  of  k.  Coupling  between  the  phonon  field 
and  the  adbond  is  brought  about  by  the  position  operator  u  of  M,  which  equals 
the  displacement  field  at  the  origin.  Explicitly, 

s  ■  l  <5Sv5rkT,i  (ak,  +  a-kS>  ?ks  <2-6) 

ks  * 


in  terms  of  the  volume  v  of  a  unit  cell  and  the  quantization  volume  V. 
Combining  everything  then  gives  for  the  Hamiltonian  of  the  adbond  plus  crystal 


H  *  H  +  H  -  RS  (2.7) 

a  p 

with  the  abbreviations 

R  «  u-e  ,  S  *  dV./dz  .  (2.8) 

-  -z  l 


It  will  appear  in  the  next  section  that  it  is  advantageous  to 
density-operator  formalism,  so  the  equation  of  motion  is 


wul  w  x  in  a 


iJ*dt  *  lH’pl  *  P+  *  P  »  Trp  «  1 


(2.9) 


From  Eq.  (2.6)  we  see  that  the  interaction  between  adbond  and  crystal,  -RS,  is 
linear  in  the  creation  and  annihilation  operators,  which  implies  that  -RS  can 
only  induce  transitions  between  crystal  states  which  differ  by  not  more  than  one 
phonon.  In  other  words,  by  retaining  only  the  linear  term  in  the  expansion 
(2.1),  we  discard  multiphonon  processes,  which  can  be  justified  as  long  as  any 
set  of  two  adbond  states  is  resonantly  coupled  by  a  single-phonon  process. 
Physically  this  implies  that  any  level  separation  |uk  -  uk  |  must  be  smaller  than 
the  cut-off  frequency  of  the  dispersion  relation,  which  is  the  Debye  frequency 

V 


III.  THERMAL  RELAXATION 

The  full  density  operator  p(t)  of  a  single  atom  and  a  large  crystal  has  not 
much  significance.  Since  we  are  interested  in  the  dynamics  of  the  adbond,  we 
consider  the  reduced  density  operator  pQ(t),  defined  as 

p(t)  -  Tr  p(t)  (3.1) 

o  p 


■r.  •  , 


where  the  trace  runs  over  all  phonon  states.  On  the  other  hand,  the  crystal  is 
regarded  as  a  large  reservoir  at  temperature  T,  for  which  the  density  matrix 
equals 


Pp  =  exp(-pHp)[Trpexp(-3Hp)]  1  (3.2) 

with  p  *  (kT)  * .  Coupling  between  the  adbond  and  the  heat  bath  then  gives  rise 
to  the  thermal  relaxation  of  p  (t),  and  eventually  p  (t)  reaches  a  steady  state 


lim  p  (t) 
t-*»  ° 


(3.3) 


in  which  the  atomic  bond  and  the  crystal  are  in  thermal  equilibrium. 

It  is  a  standard  procedure  in  relaxation  theory  to  derive  an  equation  for 

22 

the  reduced  density  operator  pQ(t).  Both  projection  techniques  and  a 

23  24 

reservoir  approach  ’  yield  the  same  result,  provided  that  equivalent 

25 

approximations  are  made.  Here  we  will  briefly  summarize  the  derivation  of 
Ref.  24  and  adopt  a  Liouville-operator  notation,  which  will  allow  a  concise 
formulation.  A  Liouvillian  L  is  related  to  a  Hamiltonian  H  according  to 

Lp  =  K'^H.p]  (3.4) 


which  defines  the  action  of  L  on  an  arbitrary  operator  p  in  Hilbert  space.  Thei 

one  regards  p  as  a  vector  in  Liouville  space  and  L  as  a  linear  operator  in  that 

26  2 
space.  Consequently,  an  N-dimensional  Hilbert  space  generates  an  N  - 

dimensional  Liouville  space.  Matrix  representations  of  Liouvillians  can  be 

constructed  in  the  very  same  way  as  for  operators  in  Hilbert  space.  A  theorem 

that  follows  immediately  from  Eq.  (3.4)  is 

exp(Ls)p  =  exp(Hs/h)pexp( -Hs/K)  (3.5) 


where,  for  instance,  s  *  it. 


Reservoir  theory  starts  with  a  transformation  of  the  equation  of  motion 
(2.9)  to  the  interaction  picture,  where  the  coupling  -RS  is  considered  as  the 
interaction.  With 


p(t)  «  exp[ i(L  +  L  )t]p(t) 
a  p 


(3.6) 


Eq.  (2.9)  can  be  written  as 


i*  dt  P(t)  "  -w(fc)p(t) 


(3.7) 


where  W(t)  is  defined  as  the  Liouvillian 


W(t)p  «  [ (exp( i ( L  +  L  )t)(RS)),p] 
a  p 


(3.8) 


for  an  arbitrary  p.  As  the  initial  condition  we  choose 


P(0)  *  p  (0)p 
o  p 


(3.9) 


Then  the  identity 


Tr  (W(t)p(O))  -  0 


(3.10) 


follows  from  the  fact  that  pp  conanutes  with  Hp,  and  from  the  explicit  form  of  R, 
Eqs.  (2.6)  and  (2.8),  which  gives  TrpRpp  *  0.  If  we  now  integrate  Eq.  (3.7), 
take  the  trace  over  the  phonon  states,  and  use  Eq.  (3.10),  we  find 


dt  PQ(t)  *  "A  2TrD  [  dT  W(t)W(t-r)p(t-x) 
“  7  n 


(3.11) 


which  is  an  exact  integral  of  the  equation  of  motion. 

Obviously,  the  integral  in  Eq.  (3.11)  is  awkward,  and  a  series  of 
approximations  has  to  be  made  in  order  to  obtain  a  manageable  expression, 
crucial  step  which  has  to  be  made  is  the  factorization 


p(t-t)  -  p  (t-t)p 
o  p 


(3.12) 


m 


in  the  integrand  of  Eq.  (3.11).  Then  we  substitute  the  explicit  form  of  W(t), 
use  Eq.  (3.5)  several  times,  and  transform  back  to  the  Schrodinger  picture  with 
the  inverse  relation  of  Eq.  (3.6).  Finally  we  obtain 


i  ^  PQ(t)  =  LaPo(t)  -  ^  Ls  J^dx  exp(iLaT)Lc(T)po(t-T) 
where  L_  and  L  (t)  are  defined  by 

J  C 

Lgp  -  tS,p} 

Lc(t)p  -  G(t)Sp  -  G( x  )*pS 


(3.13) 


(3.14) 


which  involves  the  reservoir  correlation  function 

G( x )  =  2irh"2Tr  Rp  exp(iL  x)R  .  (3.15) 

P  P  P 

Inspection  of  Eq.  (3.13)  shows  that  the  coupling  to  the  crystal  only  enters 
through  the  function  G(t).  Approximation  (3.12)  is  sufficient  to  separate 
reservoir  and  adbond  operators,  and  the  presence  of  the  solid  can  be  accounted 
for  by  a  single  function  G(x). 

Of  paramount  importance  for  the  validity  of  a  reservoir  approach  is  that 
G(t)  decays  to  zero  (for  x-**»)  sufficiently  quickly.  Due  to  the  interference 
between  the  many  phonon  modes,  the  typical  decay  time  will  be  of  the  order  of 
the  inverse  cut-off  frequency  of  the  spectrum,  which  is  for  a  crystal.  This 
will  become  more  clear  in  the  next  section.  Then  Lc(x)  deviates  from  zero  only 
for  x  —  0,  and  the  most  crude  approximation  would  be  to  replace  it  by  an 
operator  which  is  proportional  to  a  delta  function.  It  is  easy,  however,  to  do 
much  better.  The  net  effect  of  the  coupling  of  the  adbond  to  the  crystal  is  a 
damping  on  a  time  scale  which  equals  the  inverse  linewidth.  Since  this  is  very 
slow  in  comparison  with  u»_  \  we  conclude  that  in  the  interaction  picture,  p(t-x) 


varies  negligibly  on  a  time  scale  ,  and  therefore  we  can  replace  p(t-x)  by 
p(t)  in  Eq.  (3.11).  This  is  the  Markov  approximation,  and  the  equivalent  in 
the  Schrodinger  picture  reads 


PQ(t-x)  -  exp(iLaT)pQ(t)  (3.16) 

which  can  be  inserted  into  Eq.  (3.13).  Then  we  can  replace  the  upper 
integration  limit  by  infinity,  which  yields  for  the  equation  of  motion 

1  ?n(t)  *  ‘  ir)p  (t)  .  (3.17) 

dt  o  a  o 

Relaxation  of  the  adbond  due  to  the  coupling  to  the  phonon  reservoir  is  now 
conpletely  accounted  for  by  the  time* independent  Liouvillian 

T  =  L  |  dx  exp(-iL  x)L  (x)exp(iL  x)  .  (3.18) 

2ir  s  JQ  r  a  c  r  a 

An  expansion  of  T  in  matrix  elements  with  respect  to  adbond  states  will  be  given 
in  Section  V. 

If  we  insert  the  expression  for  L_  and  L  (x)  ir.  Eq.  (3.18),  and  use 

3  C 

relation  (3.5)  for  the  exponentials,  then  we  find  that  T  can  alternatively  be 
cast  in  the  form 


Tp  -  Ls(Qp  -  pQf) 

where  the  Hilbert-space  operator  Q  is  defined  as 
Q  *  J"dx  G(x)exp(-iLax)S  . 


(3.19) 


(3.20) 


The  advantage  is  that  Q  only  involves  a  single  exponential.  Furthermore,  we 
directly  find  from  the  representation  (3.19) 


which  guarantees  the  conservation  of  trace  and  Hermiticity  in  the  time  evolution 

of  o  (t). 
o 

IV.  RESERVOIR  CORRELATION  FUNCTION 

Properties  of  the  crystal  which  affect  the  time  evolution  of  the  adbond  are 
*11  embodied  in  the  correlation  function  G(t).  In  order  to  illustrate 
quantitatively  the  behavior  of  G(t),  we  adopt  a  simple  Debye  model,  which  will 
already  reveal  the  most  salient  features  of  a  crystal  correlation  function.  In 
this  model,  the  dispersion  relation  is  taken  to  be 


w(k)  *  ckH(ujp-ck) 


(4.1) 


with  c  the  speed  of  sound  and  H  the  unit-step  function.  The  expression  (3.15) 
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is  easily  evaluated  with  standard  techniques.  It  is  elucidating  to  subdivide 
G(t)  in  a  spontaneous  and  a  stimulated  part,  according  to 


G(t)  *  G(t)  +  G(t) 

sp  st 


(4.2) 


where  G(t)  is  by  definition  G(t)  for  T  »  0.  Hence  G(t)  is  independent  of  T 
sp  sp 

and  accounts  for  relaxation  of  the  adbond  due  to  the  presence  of  the  crystal, 
irrespective  of  its  temperature.  We  find  for  this  spontaneous  part 


G(t) 


sp 


j  (1  +  iuijjT  )exp(  -  iujpT )  -  1 

Cud  - : — 72 - 

(v} 


(4.3) 


with  the  parameter  £  given  by 


C  -  3*/*M 


(4.4) 


The  real  and  imaginary  part  of  G(t)  are  plotted  in  Fig.  1.  For  the  stimulated 

sp 


contribution  we  obtain 
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G(t  ) 


st 


dw  n(w)wcos(u)T ) 


(4.5) 


which  is  plotted  in  Fig.  2  for  two  values  of  the  temperature.  Here  n(u»)  is  the 
average  number  of  phonons  in  mode  u>,  which  is  explicitly 


n(w)  -  [exp($rtw)  -  1] 


(4.6) 


Figs.  1  and  2  illustrate  that  G(t)  tends  to  zero  indeed  on  a  time  scale  i»q~  , 
although  the  decay  is  not  exponential. 

From  Eq.  (3.20)  it  follows  that  the  relaxation  constants  (matrix  elements 
of  ")  are  determined  by  the  Fourier-Laplace  transform  of  G(t) 


S(u>) 


dr  exp(  iwt)G(-r) 


(4.7) 


with  w  equal  to  a  transition  frequency  uk 
readily  derive 


u>. . 

J 


From  Eqs.  (4.3)  and  (4.5)  we 


S(w)sp  ■  Cuijj  \jH(w)H(uip-w)  -  ix  3(<Up  +  wlog|  l-Wp/w|  )  (4.8) 

C(w)st  *  Cwp  3|w|n(  ju)|  )H(wp  -  M) 

+  iir'1^'3?  r°  du'  n(u)'  )  (4.9) 

'0  w  -  a)' 


where  P  stands  for  principal  value.  These  functions  are  plotted  in  Figs.  3  and 
4.  For  negative  values  of  <d  the  real  part  of  6(w)gp  vanishes  identically,  which 
implies  that  this  contribution  to  T  represents  a  decaying  part.  The  thermal 
part  obeys  the  relation 


0(  -w) 


st 


(4.10) 


V.vIvV 


which  indicates  that  upward  and  downward  transitions  between  two  levels,  due  to 
the  finite  temperature  of  the  crystal,  have  the  same  rate  constant.  Hence  the 
genuine  thermal- equilibrium  distribution  is  brought  about  by  the  spontaneous 
part  of  the  relaxation  operator,  which  is  temperature  independent. 

Setting  u>  equal  to  zero  yields  immediately 

t)  *  Culp  ^(-i  +  irkT/Jtap)  (A.  11) 

which  can  be  regarded  as  a  measure  of  the  strength  of  G(t).  For  kT>Kui^  the 
stimulated  transitions  tend  to  dominate  the  spontaneous  decay. 


V.  MATRIX  ELEMENTS 

Solving  the  equation  of  motion  (3.17)  for  a  particular  configuration  of 
adbond  states  requires  an  expansion  in  matrix  elements.  Taking  the  k,l-th 
matrix  element  gives 

^  <k|po(t)|t>  -  -i<k|(LaPo(t))|*>  -  <k|(rpo(t))|l>  .  (5.1) 

It  remains  to  express  the  right-hand  side  in  matrix  elements  of  pQ(t).  From  the 
definition  (3. A)  of  L  and  the  expression  (2. A)  of  H  ,  we  easily  find  for  the 
first  term 

<k|(Lap)|t>  ■  Akl<kjpjt>  (5.2) 

for  an  arbitrary  p,  where 

Set  '  “k  •  “t  (5-3) 

is  the  level  separation  between  the  adbond  states  |k>  and  |t>. 

Slightly  more  complicated  is  the  evaluation  of  the  relaxation  operator  T. 
First  we  notice  that  the  exponential  in  Eq.  (3.20)  can  be  expanded  as 


which  yields  the  representation  for  Q 


o-il  5<atk)pkspt  (5-5) 

kt 

in  terms  of  the  Fourier-Laplace  transform  of  the  correlation  function.  Then  we 
insert  Eq.  (5.5)  into  Eq.  (3.19)  and  use  the  closure  relation 

£  ?i  «  1  (5.6) 

i 

which  finally  gives 

<k|(rp)|t>  «  i  £  <ckmmn<n|p|a>  +  ctmmn<klpln>) 

mn 

-  *  l  <cnlk»  +  (5-7) 

mn 

in  terms  of  the  parameters 

c,  „  *  <k | S 1 t><m|s|n>G( A  )  .  (5.8) 

Equation  (5.7)  expresses  the  relaxation  operator  T  in  matrix  elements  of  the 
derivative  of  the  potential,  S  *  dV^/dz,  and  G(w),  evaluated  at  the  frequencies 
The  combination  of  Eqs.  (5.1),  (5.2)  and  (5.7)  turns  the  equation  of 
motion  into  a  simple  set  of  linear  first-order  differential  equations,  which  can 
be  solved  directly  once  a  potential  and  an  initial  state  pQ(0)  are  specified. 

VI.  ABSORPTION  SPECTRUM 

A  density  operator  is  not  directly  amenable  to  observation  in  an 
experiment.  One  way  of  measuring  properties  of  an  adsorbed  atom  is  by  probing 


the  system  with  a  weak  monochromatic  laser.  A  field  with  intensity  1^  (energy 
per  unit  of  time  that  passes  through  a  unit  area,  perpendicular  to  the  direction 
of  propagation)  and  polarization  is  scanned  over  the  resonances  of  the 
system,  and  the  power  absorption  I(uj^)  (absorbed  energy  per  unit  of  time)  is 
measured  as  a  function  of  the  laser  frequency  Photons  are  absorbed  from  the 

field  by  the  joint  system  of  atom,  crystal  and  their  interaction  (although  the 
atom  and  the  crystal  separately  are  both  assumed  to  be  transparent),  so  that  the 
absorption  profile  exhibits  the  details  of  the  coupling  between  the  adbond  and 
the  phonon  field,  rather  than  the  properties  of  the  adsorbed  atom  (the 
potential)  alone. 

The  interaction  between  the  adsorbate  and  the  radiation  field  is 
established  by  a  dipole  coupling.  Since  the  motion  of  the  atom  is  restricted  to 
the  z-direction,  the  dipole-moment  operator  p  can  be  written  as  pez,  with  p  an 
operator  in  the  Hilbert  space  spanned  by  the  eigenstates  |k>  of  the  adbond 
Hamiltonian  Ha«  With  E  the  electric  component  of  the  radiation  field,  we  can 
include  the  interaction  in  the  equation  of  motion  (2.9)  with  the  substitution 

H  -  H  -  pe  -E  .  (6.1) 

-z  - 

The  probe  beam  is  considered  to  be  very  weak,  so  that  we  can  compute  the  power 
absorption  with  the  Golden  Rule.  From  the  Appendix  of  Ref.  1  we  copy  the  formal 
result 

Ko^)  ■  ILuL(e0Kc')  1|ez-cLJ2Re  j"dT  expdta^-OTrptpf-O.p]  (6.2) 

where  c'  is  the  speed  of  light  and  p  is  the  thermal-equilibrium  density  operator 
of  the  entire  system,  but  without  the  interaction  -p-E.  The  T-dependence  of 


Vi(")  Is  given  by  y(x)  -  exp(iLt)p.  We  remark  that  the  commutator  gives  rise  to 
two  terms,  where  p(x)p  (pp(t))  represents  stimulated  photon  absorption 
(emission).  The  net  absorption  is  the  balance  between  the  loss  and  gain  term 
for  photons  in  the  laser  mode. 

Evaluation  of  expression  (6.2)  starts  with  the  identity 


Trp[ii(t) ,p)  »  Tr pexp( - iLr ) [ y , p ]  . 


(6.3) 


Then  we  introduce  the  quantity 


D(t)  *  exp(-iLx)[M,p] 


(6.4) 


which  is  an  operator  in  the  entire  phonon  and  adbond  Hilbert  space.  With  Dq(t) 


TrpD(t),  the  adbond  part  of  D(x),  we  can  write  Eq.  (6.3)  as 


Trp[ p( r ) , p]  =  Tr  pD  (t)  . 

a  o 


(6.5) 


Then  we  notice  that  D(t)  obeys  the  differential  equation 


.d 


i^-D(x)  =  LD(t) 


(6.6) 


with  initial  condition 


D(0)  *  (p,p  ]p  (6.7) 

o  p 

where  we  have  assumed  that  p  -  p  p  .  The  equation  of  motion  (6.6)  is  identical 

-op 

to  Eq.  (2.9),  which  is  sometimes  referred  to  as  the  quantum- regress ion 
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theorem.  This  implies  that  the  time  regression  of  the  dipole  correlation 
function  D(t)  is  identical  to  the  time  evolution  of  the  density  operator  p(t). 
This,  in  turn,  is  tantamount  to  the  statement  that  the  dynamical  properties  of 
the  system  are  reflected  in  the  frequency  dependence  of  the  absorption  profile, 
which  elucidates  the  significance  of  a  measurement  of  I(w,  ). 


From  Eq.  (6.5)  we  see  that  we  only  have  to  solve  Eq.  (6.6)  for  the  reduced 
correlation  function  Dq(t),  which  can  be  done  along  the  very  same  lines  that  led 
to  Eq.  (3.17).  Hence  we  can  immediately  write  down  the  solution,  which  is 

D  (t)  *  exp(-i(L  -iDOlp.p  1  •  (6.8) 

o  a  o 

Substitution  into  Eq.  (6.2)  and  performing  the  t  integration  then  yields 

I(wL)  *  VL(eo*c,) "1I2«'£LI2r*tV  —  f [ M , PQ ]  (6.9) 

in  terms  of  the  steady- state  solution  pq  from  Eq.  (3.17),  which  obeys 

(L  -  iDp  «  0  ,  Tr  p  -  1  ,  5  +  -  p  .  (6.10) 

a  o  a  o  o  o 

From  conservation  of  trace  in  the  time  evolution  of  p  (t),  it  follows  that 

o 

lim  D  (t)  »  p  Tr  [y,P  ]  =  0  (6.11) 

o  o  a  o 

T-+* 

which  proves  that  the  principal  part  in  the  upper  integration  limit  r  *  • 
vanishes.  Equation  (6.10)  is  easily  solved  for  any  level  configuration.  Then 
we  insert  pQ  into  Eq.  (6.9),  which  determines  the  absorption  profile  in  terms  of 
a  matrix  inversion. 

VII.  LINE  SHAPE 

Expression  (6.9)  represents  the  complete  absorption  profile  as  a  function 
of  In  order  to  disentangle  the  contributions  from  the  various  adbond 

resonances,  we  consider  the  situation  of  two  levels  |l>  and  |2>  with  separation 
u>2  ~  >  0,  and  we  scan  the  laser  over  this  resonance.  The  resulting 

profile  Kwj^)  is  then  termed  the  line  shape  of  this  particular  transition.  With 
the  expansion  in  matrix  elements  from  Section  V,  we  can  write  down  immediately 
the  equation  of  motion  as  a  linear  set  of  differential  equations.  With  p^  * 

<1 1  po J 2>  and  so  forth  for  the  other  p^,  we  btain 


1  dt  p22 


_ia21p22  +  U12pll 


(7.1) 


1  dt  P11 


ia21p22  ‘  ia12pl 1 


(7.2) 


.  d_ 

1  dt  P21 


=  (wo-in)p21  +  in  exp(2i(pg)p12 


(7.3) 


,  d 

1  dt  p12 


inexp(-2i<»s)p21  -  (uo+in  )pj 


(7.4) 


where  we  introduced  the  abbreviations  for  the  relaxation  parameters 


akt  -  |  <k | S  |  i>|  ReS(Aj^) 


(7.5) 


n  *  i|<2|s| l>|2(C(u)o)  +  5(-(uq)  ) 


(7.6) 


and  neglected  the  couplings  due  to  the  small  diagonal  matrix  elements  <l|S|l> 
and  <2|S|2>.  Furthermore,  the  phase  ijig  of  <2 1 S 1 1  >  is  defined  by 


< 2 1 S j 1 >  *  |<2|s|l>|exp(i*_) 


(7.7) 


Of  course,  we  can  choose  the  phase  between  |l>  and  |2>  in  such  a  way  that  <J>g  * 
0,  but  it  will  turn  out  to  be  advantageous  to  postpone  the  definition  of  the 
relative  phase. 

The  steady-state  solution  of  the  set  (7.1)  -  (7.4)  is  readily  found  to  be 


"2  *  a21/U21  +  *12)  •  nl  '  *12/(*21  *  a12> 


(7.8) 


for  the  populations  n^  ■  p^,  whereas  the  coherences  p^2  and  p2^  vanish.  From 
the  set  (7.1)  -  (7.4)  we  read  off  the  matrix  representation 


WL  +  U12 


-  L  +  ir 
a 


“l  '  wo  A  in 

-inexp(-2i^  ) 

S 


-in  exp(2i*  ) 
s 

“l  +  wo  +  in* 


on  the  basis  |2><2|,  |l><lj,  |2Xl|,  |l><2|.  The  inverse  of  this  matrix  is  then 
substituted  into  Eq.  (6.9).  Matrix  elements  of  the  dipole-moment  operator  are 
not  necessarily  real.  We  write 


<2|p|l>  *  | <2 1 y |  l>jexp(iijO  (7.10) 

and  now  we  choose  the  phase  of  the  wave  function  |l>  in  such  a  way  that  t|>  *  $  . 

M  S 

Combining  everything  then  gives  for  the  absorption  line 


K“L>  -  :lB(51  • 


4ui  at.  Re 
oh  n 


(a)2  -  uj?  +  2u>  Imq)2  + 

O  L  O 


4u2(Ren)2 


(7.11) 


where  the  Einstein  coefficient  B  for  stimulated  transitions  between  |l>  and  |2> 
is  defined  as 


B  "  )  1|<2lpl1>?z‘-L^2  • 


(7.12) 


Due  to  the  fact  the  Ren  and  Imr|  are  not  necessarily  small  in  comparison  with  idq, 
the  line  is  not  a  simple  Lorentzian.  A  few  examples  of  l(^)  are  drawn  in  Figs. 
S  and  6. 


VIII.  SECULAR  APPROXIMATION 


A  great  simplification  arises  if  the  relaxation  constants  c^^^  are  small 
in  comparison  with  the  transition  frequencies  A^  of  the  free  evolution  of  the 
adbond.  In  spectral  terms  this  means  that  the  widths  and  shifts  of  lines  are 
small  in  comparison  with  their  central  frequency,  so  ) n  | <<ot>o  in  the  notation  of 
the  previous  section.  Then  the  Liouvillian  will  dominate  the  time  evolution 
of  p  (t)  in  Eq.  (3.17),  and  the  coupling  between  eigenvectors  of  L  with 
different  eigenvalues  can  be  neglected.  Every  frequency  A^  gives  rise  to  a 


single  line,  and  the  time  regression  of  the  correlation  function  does  not  couple 
anymore  between  different  lines.  In  this  approximation,  every  line  evolves  in  a 


secular  fashion,  which  explains  the  origin  of  the  name  "secular"  for  this 
Uniting  approximation. 
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Mathematically  this  implies  that  we  can  decouple  the  time  evolution  of  the 

coherences  <k|p  (t)|i>,  k  *  l,  from  the  equation  for  the  populations  n,  (t)  = 

<k|p  ( t ) | k> ,  since  the  free  evolution  of  <k I p  (t)|i>  is  proportional  to  exp(- 
■  o  1  1  o  1 

^t),  whereas  n^(t)  evolves  with  eigenvalue  zero  (A^  M  0).  From  the  matrix 
representation  of  section  V  it  then  follows  that  the  coherences  decay 
exponentially  to  zero,  and  that  the  time  evolution  of  the  populations  is 


governed  by  the  master  equation 


dT  "k  ’  £(n 


„  3.  .  ,  “  n.  d.  . ) 

l  ik  k  ki 


(8.1) 


where  only  the  relaxation  constants  a^  from  Eq.  (7.5)  appear.  Equation  (8.1) 
is  a  simple  gain- loss  balance  for  the  population  of  level  |k>,  and  we  can 
interpret  the  term  n^a^  as  the  rate  of  transitions  from  | t>  to  |k>,  due  to 
single-phonon  emission  into  the  crystal  (w^  >  ui^),  or  energy  absorption  from  the 
solid  (w^  <  ui^)-  The  set  of  differential  equations  (8.1)  should  be  accompanied 
by  the  constraint 


(8.2) 


which  expresses  conservation  of  trace. 

From  Eq.  (8.1)  it  follows  that  the  decay  constant  for  level  |k>  can  be 
written  as 


S-I 


(8.3) 


which  turns  the  master  equation  into 


We  notice  that  in  the  secular  limit,  the  relaxation  constants,  which  determine 


the  level  populations,  are  real.  The  coherences,  however,  still  acquire 
contributions  from  the  imaginary  parts  of  the  cic^mn's*  which  amounts  to  the  line 
shifts.  An  advantage  of  the  secular  limit  is  that  the  imaginary  parts  can  be 
accounted  for  by  attributing  an  effective  level  shift  to  each  frequency 
according  to 


“k  =  “k 


+  J^Imc 


ktfck 


(8.5) 


Then  we  consider  this  transformation  done,  and  suppress  the  tilde  henceforth. 

For  further  use,  we  remark  that  in  the  secular  limit  the  relaxation  operator  can 
be  expressed  entirely  in  projection  operators.  We  obtain 

rp  -  *  l  akt(V  +  »pk  -  lWk<>>  (8-6) 

kl 

as  can  be  checked  by  inspection. 

Obviously,  the  secular  approximation  is  not  suitable  for  the  evaluation  of 
the  details  of  a  line  shape,  since  it  turns  every  line  into  a  Lorentzian.  In 
the  remainder  of  this  paper  we  shall  consider  the  effect  of  irradiation  with  a 
strong  laser.  Then  the  dominant  features  of  this  system  will  be  determined  by 
the  driving  field,  rather  than  by  the  minor  details  of  the  coupling  to  the 
phonon  field.  Hence  we  can  adopt  the  representation  (8.6)  for  the  relaxation 
operator  from  here  on. 


IX.  COHERENT  EXCITATION 

Probing  the  adsorbate  with  a  weak  laser  does  not  alter  the  dynamics  of  the 
system.  In  a  different  application  of  lasers  in  these  configurations,  one 


>W\  . 

A  A  vV-V 


deliberately  tries  to  modify  or  affect  the  behavior  of  the  system  by  driving  a 
particular  transition  with  a  strong,  resonant  field.  The  coherent  nature  of  a 
laser  field  (well-defined  phase  in  its  time  evolution)  provides  essentially  a 
different  excitation  mechanism  than  a  coupling  to  a  thermal  bath,  like  the 
phonon  reservoir.  Along  with  the  fact  that  the  laser  power  can  be  very  high, 
this  then  opens  the  possibility  to  drive  the  system  away  from  thermal 
equilibrium.  The  presence  of  the  radiation  will  tend  to  maintain  a  certain 
distribution  of  the  population  over  the  different  levels,  which  has  to  compete 
with  the  thermal  relaxation.  Therefore,  it  can  be  anticipated  that  it  should  be 
feasible  to  actually  change  the  dynamics  of  the  adbond  if  the  laser  is 
sufficiently  intense.  As  mentioned  in  the  Introduction,  the  first  goal  has  been 
to  enhance  the  desorption  by  resonant  excitation  of  the  adbond. 

Radiative  transitions  between  vibrational  adsorbate  states  have  been 

28*32 

studied  extensively,  both  in  the  weak-field  case  and  for  strong 

fields.**’  ^  ^  In  the  present  paper  we  summarize  and  extend  our  own 
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approach,  which  is  valid  for  arbitrary  laser  power  and  includes  the  effect 

of  the  laser  linewidth.  We  shall  again  work  with  a  Liouville  notation,  which 
allows  a  concise  formulation.  Suppose  that  the  laser  frequency  is  in  close 
resonance  with  a  single  transition  of  the  adbond  only,  and  that  it  couples  a 
ground  state  |g>  (not  necessarily  the  lowest  state  of  the  adbond)  with  an 
excited  state  ]e>,  which  are  separated  by 


b> 


u>  -  u>  >0 

e  g 


(9.1) 


The  idea  is  that  |e>  is  one  of  the  high- lying  states  in  the  potential  well, 
which  is  very  unlikely  to  be  populated  by  thermal  excitation.  Then  this 
configuration  automatically  excludes  resonant  coupling  between  other  levels. 


Radiative  transitions  only  occur  between  |e>  and  |g>,  but  the  thermal  coupling 
remains  present  between  all  levels. 

The  electric  component  of  the  laser  field  at  the  position  of  the  adbond  is 
given  by 


E(t)  *  Eq  Re  e^exp { - i [ +  $(t)]} 


(9.2) 


where  $(t)  is  a  real-valued  stochastic  process,  which  is  responsible  for  the 

3( 

broadening  of  the  laser  line  around  its  central  frequency  It  can  be  shown 

that  the  response  of  the  system  is  quite  insensitive  for  the  stochastic  details 
of  $(t),  provided  that  we  restrict  the  description  to  single-mode  fields.  We 
shall  take  $(t)  as  the  independent -increment  process, ^  which  covers  the  more 
familiar  Gaussian  white-noise  and  random-jump  processes  as  special  cases.  The 
multiplicative  stochastic  differential  equation  for  this  diffusion  process  has 
been  solved  in  the  Appendix  of  Ref.  41,  where  it  turns  out  that  the  spectral 
profile  of  the  laser  is  a  Lorentzian,  whose  half -width  at  half-maximum  is 
denoted  by  X. 

The  dipole  coupling  between  the  adbond  and  the  laser  field,  Eq.  (6.1),  now 
attains  the  explicit  form 

-U-E(t)  *  -ihn|e><g|exp{-i(w^t  +  $(t)J)  +  Hermitian  conjugate  (9.3) 


where  the  usual  rotating-wave  approximation  has  been  made,  and  permanent 
dipole -moments  have  been  omitted.  The  strength  of  the  coupling  is  determined  by 
the  Rabi  frequency 


0  *  *  1Eo|<e|u-eL|g>|  • 


(9.4) 


Then  the  equation  of  motion  for  'he  reduced  density  operator  becomes 
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i#-7TP  (t)  ■  tH  -  p-E(t),p  (t)]-i#rp  (t) 

at  o  a  -  -  o  o 


(9.5) 


which  now  contains  the  interaction  with  the  external  field  explicitly. 


a  convenient  way  to  eliminate  the  fast  oscillations  with  the  optical 


frequency  from  the  Hamiltonian  is  by  the  introduction  of  the  stochastic 
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unitary  transformation 


o(t)  *  exp{-i[w^t  +  $(t)]L  }pQ(t) 


(9.6) 


which  contains  the  Liouvillian 


L  p  *  [P  ,p]  • 

g  g 


(9.7) 


It  is  easy  to  check  that  the  transformation  only  affects  the  coherences,  so  that 


we  have 


<k|  o( t )  |k>  *  n^(t) 


(9.8) 


for  every  level  |k>.  With  some  algebra,  we  find  the  transformed  equation  of 


motion  to  be 


ijr  o(t)  *  (L.  +  ♦(t)L  -  ir)o(t) 

dt  d  g 


(9.9) 


with  LdP  *  K  [Hd,p]  and 


Hd  -  Ha  +  tajP  -  i#n(  |e><g|  +  |g><ej) 


(9.10) 


the  dressed-adbond  Hamiltonian.  This  Hd  has  the  significance  of  representing 


the  free  adbond  (H  ),  the  free  evolution  of  the  single  laser-mode  and  their 
a 


dipole  coupling.  Phonon  transitions  (T)  and  phase  fluctuations  ($(t)L  )  couple 

© 


the  eigenstates  of  Hd»  which  both  give  rise  to  relaxation  of  o(t). 


- V 


The  appearance  of  the  time  derivative  $(t)  of  the  stochastically 
fluctuating  phase  turns  the  equation  of  motion  for  o(t)  into  a  multiplicative 
stochastic  differential  equation,  and  the  density  operator  o(t)  into  a 
stochastic  process.  Only  the  average  over  many  realizations  of  the  process  t(t) 
can  have  relevance,  and  we  write 

n(t)  »  |o(t)}  (9.1D 

where  the  brackets  indicate  the  average.  Then  Eq.  (9.9)  is  easily 

solved  for  its  average,  and  we  obtain  the  equation  of  motion  for  n(t) 

i~n(t)  *  (Ld  -  iw  -  ir)n(t)  (9.12) 

where  the  effective  relaxation  operator,  which  accounts  for  the  finite  laser 
linewidth,  is  given  by 

W  =  AL2  .  (9.13) 

8 

This  operator  can  alternatively  be  represented  as 

Wp  -  A(P  p  +  pP  -  2P  Tr  P  p)  (9.14) 

g  g  gag 

which  is  reminiscent  of  the  structure  of  T,  Eq.  (8.6).  In  Eq.  (9.12)  we  can 

incorporate  the  effect  of  the  laser  linewidth  by  the  substitution  a  -*■  a  +  2\ 

gg  gg 

in  the  definition  of  T,  although  in  general  this  is  not  the  correct  procedure, 
as  we  will  see  in  due  course. 

X.  DRESSED  STATES 

Both  the  laser  linewidth  and  the  coupling  to  the  phonon  reservoir  give  rise 
to  a  damping  of  the  free  evolution  of  the  dressed  adbond,  which  is  represented 
by  the  Liouvillian  in  Eq.  (9.12).  This  is  the  analogue  of  from  Eq. 
(3.17),  pertaining  to  a  field- free  system.  In  order  to  illuminate  the  physical 
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picture,  we  diagonalize  the  Hamiltonian  H, .  First  we  rewrite  Eq.  (9.10)  as 


H ,  *  )  rtw.P.  +  J)Hu)  +  w  +  uj.  )(P  +  P  ) 
d  i  i  i  egLeg 


i*e,g 


-  i*A(Pe  -  Pg)  -  ilta(|eXg|  +  |g><e|) 


(10.1) 


with  A  *  -  a»o  the  detuning  from  resonance.  Diagonalization  of  is  now 


trivial,  and  we  find  the  eigenvalue  equations  to  be 


Hd|k>  «  llaijJlO  ,  k  *  e,g 


(10.2) 


Hd|i>  *  +  00  +  *F  £2  *  )  j  ±> 


(10.3) 


in  terms  of  the  generalized  Rabi  frequency 


n'  =  a(  l  +  o2/aV 


(10.4) 


Hence  the  dressed  states  Jk>  with  k  *  e,g  are  equal  to  the  adbond  states  |k>, 
and  they  have  the  same  eigenvalue,  whereas  the  states  |e>  and  |g>  form  the 


linear  combination 


|+>  =  jg>sini©  +  |e>cosJ© 


(10.5) 


|->  *  |g>cosi©  -  |e>sin|0 


(10.6) 


which  is  parametrized  with  the  angle 


9  ■  arctan  (Q/A) . 


(10.7) 


The  position  of  the  dressed  states  with  respect  to  the  adbond  states  is 


illustrated  in  Fig.  7. 


XI .  STEADY  STATE 
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Solving  Eq.  (9.12)  for  a  transient  regime  [0,t]  has  little  significance, 
since  a  preparation  of  a  specific  initial  state  11(0)  is  practically  not 
feasible.  Any  state  n(0),  however,  will  relax  to  the  same  long-time  solution  IT 


on  a  time  scale  of  the  order  of  r 
the  equation 


(L.  -  iW  -  ir)H  «  0 


-1 


Therefore,  we  consider  in  some  more  detail 


(11.1) 


Due  to  the  presence  of  the  driving  field,  the  equations  for  the  populations  11^  ■ 
<k|n|t>  will  not  necessarily  decouple  anymore  from  the  equations  for  the 
coherences  between  different  states.  In  general,  all  matrix  elements  ofTT(t) 
couple,  but  it  turns  out  that  in  the  steady  state  the  coherences  for  k  *  i. 
vanish,  provided  that  (k,t)  *  (e,g).  Taking  matrix  elements  of  Eq.  (11.1)  shows 
that  the  relation  between  the  populations  only  involves  the  coherence  between 
the  driven  transition  |e>  -  |g>,  as  could  be  expected.  Using  TT+  *  TI  gives  the 
important  relation  between  the  real  and  imaginary  parts 


(Ag  +  Ag  +  2X)Re<e|IIjg>  *  -2AIm<e JTT J g> 


which  enables  one  to  combine  the  equations  for  H  and  II  as 

g  e 


(11.2) 


a  (H  -  II  )  *  nim<e|TI|g> 

IT  © 


(11.3) 
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in  terms  of  the  parameter 


a  *  JO 
r 


5  i<A  +  A„ )  +  X 
2 _ ® _ 8 _ 

(J(A  +  A  )  +  X}2  +  A2 

e  g 


(11.4) 


With  the  aid  of  Eqs.  (11.2)  and  (11.3)  we  can  eliminate  the  coherence,  which 
finally  yields  the  set  of  equations  for  the  populations 


l  Vo.  ■  Vk  •  k  * 


(11.5) 
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\  Ha  *  a  n  -  a(n  -  n) 

l  l  le  ee  rg  e 


(11.6) 


y  Il.a.  =  A  IT  +  a  (IT  -  IT) 
L  1  8  8  rg  e 


(11.7) 


This  set  can  be  considered  as  a  steady-state  analogue  of  the  master  equation 


(8.1), 


Inspection  of  Eqs.  (11.5)  -  (11.7)  shows  that  we  can  write  the  set 


alternatively  as 


l  Vik  *  l  Vit 


(11.8) 


where  the  primed  constants  are  defined  as 


.  k,t  *  e,g  or  g,e 


akl  +  .  k,*.  »  e,g  or  g,e 


(11.9) 


In  the  same  fashion  as  when  we  identified  the  quantity  H^a^  as  the  rate  of 

transitions  from  |k>  to  | t>  due  to  single-phonon  processes,  we  can  now  interpret 

II  a  and  Ha  as  the  rates  of  stimulated  radiative  transitions  from  |e>  to  |g> 
e  r  g  r  1,0 

and  from  |g>  to  |e>,  respectively.  We  shall  see  in  Section  XIV  that  the  optical 


transitions  acquire  contributions  from  both  single-photon  and  multiphoton 


processes.  Notice  that  the  three  optical  parameters  £3,  A  and  X  only  enter  the 


equations  for  the  populations  through  the  combination  (11.4)  in  a 


Furthermore,  the  rate  constant  ar  is  linear  in  the  laser  power  (-  u),  which 


implies  that  the  number  of  radiative  transitions  increases  indefinitely  with  an 


increasing  incident  intensity.  Obviously,  the  net  absorbed  power  by  a  single 


adsorbed  atom  should  reach  a  saturation  value  in  the  limit  of  high  irradiances. 


This  will  be  shown  in  the  next  section. 


XII.  LASER  HEATING 


Radiative  excitations  of  the  adbond  occur  at  a  rate  a  I!  ,  whereas 

r  g* 


stimulated  emissions  of  photons  in  the  laser  field,  accompanied  with  )e>  -*•  )g> 
transitions,  happen  a^I^  number  of  times  per  unit  of  time.  Balancing  the  rates 
gives  an  effective  rate  of  a^OI^  "  Ee)  f°r  the  absorption.  Every  transition 
corresponds  to  an  effective  absorption  of  a  laser  photon,  so  that  the  power 
absorption  should  roughly  be  equal  to  fou^a^Ug  -  TI^).  Since  the  accumulated 
energy  in  the  vibrational  bond  must  be  constant  in  the  steady  state,  the 
absorbed  power  from  the  laser  beam  equals  the  power  flow  into  the  crystal.  This 
process  of  laser  heating  of  the  substrate  is  entirely  mediated  by  adsorbates, 
because  the  crystal  itself  was  assumed  to  be  transparent. 

In  Section  VI  we  evaluated  the  power  absorption  from  a  weak  monochromatic 
incident  field,  by  applying  the  Golden  Rule.  Obviously,  this  approach  fails  for 
strong,  finite  linewidth  radiation,  so  that  we  must  find  another  way  to 


calculate  the  power  absorption.  This  is  accomplished  by  considering  the  work 

42 


done  on  the  dipole  moment  by  the  external  field,  which  is  formally  given  by 


P(t)  -  E(t)-^<u(t)> 


(12.1) 


Here,  u(t)  -  exp(iLt)u  and  the  angle  brackets  denote  the  quantum  expectation 
value.  Transformation  to  the  Schrodinger  picture  gives 


dp. 


P(t)  -  Tr^lEC  t)  - 


(12.2) 


which  only  involves  the  reduced  density  operator  PQ(t)  of  the  adbond,  rather 


than  the  Liouvillian  L  in  Eq.  (12.1),  which  pertains  to  the  entire  system.  With 
transformation  (9.6)  we  go  to  the  o-representation,  which  yields 


P(t)  ■  ^dt  Re<el°l8>  +  +  ♦  (t))Ini<e}o}g>} 


(12.3) 
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Application  of  the  equation  of  motion  (9.9)  for  o(t)  results  in 


P(t)  *  -  itffl(Ae  +  Ag)Re<e|o|g>  +  hflwoIm<e|o|g> 


(12. A) 


which  elminates  both  the  time  derivative  and  $(t)  from  expression  (12.3). 


Equation  (12. A)  clearly  exhibits  that  power  absorption  from  an  external  field  is 


reflected  in  the  appearance  of  coherences  in  the  system. 


The  quantity  P(t)  depends  stochastically  on  time,  due  to  the  fluctuations 


in  the  laser  phase.  If  we  define  the  steady- state  power  absorption  by 


Kt^)  =  lim  {P(t)f 
t 


(12.5) 


then  we  obtain  X ( )  from  Eq.  (12. A)  with  the  substitution  o  -+  II.  Subsequently 


we  use  Eqs.  (11.2)  and  (11.3),  which  amounts  to 


2A(ul  -  on  ) 

I(V  -  *<“L  '  A  '  V 

eg  8 


(12.6) 


For  a  monochromatic  laser  this  reduces  to  I(u^)  *  ttu^a^dl  -  TI^ ) ,  as  anticipated 


in  the  beginning  of  this  section.  With  increasing  laser  linewidth  A,  the  factor 


in  curly  brackets  tends  to  uiq,  which  reflects  that  for  A  large,  the  frequency 
loses  its  significance.  Then  photons  are  considered  to  be  absorbed  in  an  |e>  ■+ 


|g>  transition,  which  correponds  to  an  excitation  of  the  sytem  with  energy  ltaQ. 


Furthermore,  we  remark  that  the  effect  of  the  laser  linewidth  in  Eq.  (12.6) 


cannot  be  incorporated  by  the  simple  substitution  a^  -*•  a^  +  2A,  as  was  the 


case  for  the  equation  of  motion  for  n(t). 


With  Eq.  (11.6)  we  can  cast  expression  (12.6)  in  the  form 


Koo,  ) 


KVe  '  l  Vu> 


(12.7) 


where  is  an  abbreviation  for  the  term  in  curly  brackets  in  Eq.  (12.6). 
Representation  (12.7)  of  the  power  absorption  only  involves  phonon-relaxation 
constants,  and  not  a^  anymore.  From  the  restriction  0  <  <  1 ,  which  holds  for 

every  population,  we  immediately  deduce  the  upper  limit 

Ku^)  <  *w£Ae  (12.8) 

for  the  laser  heating  of  the  crystal.  Although  the  rate  constant  a^  for 
stimulated  transitions  can  become  arbitrarily  large,  the  net  power  absorption 
exhibits  a  saturation,  where  the  upper  limit  is  set  by  the  phonon  rate  constants 
rather  than  by  optical  parameters.  This  can  be  understood  frcm  the  fact  that  an 
absorbed  photon  can  only  be  converted  into  thermal  energy  through  a  phonon 
transition. 

XIII.  TRANSITIONS  BETWEEN  DRESSED  STATES 

In  Section  VIII  we  adopted  the  secular  approximation  for  the  phonon- 

transition  operator  r  with  the  argument  that  the  time  evolution  of  PQ(t)  is 

dominated  by  the  f ree-evolution  Liouvillian  L’a  of  the  adbond.  This  gave  rise  to 

the  identification  of  the  rate  constants  a^  for  the  transitions  |k>  ■+  |t>, 

which  occur  at  a  rate  n^Oa^  .  For  the  strongly-driven  adbond,  however,  the 

time  evolution  of  the  system  is  governed  by  Eq.  (9.12)  for  n(t),  which  contains 

the  free  evolution  of  the  dressed  adbond  as  the  dominant  part.  Hence  it 

would  seem  to  be  more  appropriate  to  define  the  secular  approximation  with 

respect  to  the  dressed  states,  rather  than  with  respect  to  the  field-free 

states.  Essentially,  we  should  start  with  Eq.  (2.9),  include  the  interaction 

-p-E(t)  with  the  laser  field  in  the  Hamiltonian,  transform  to  the  o-picture,  and 

incorporate  the  phonon  coupling.  This  procedure  gives  an  expression  for  T  in 

the  presence  of  a  driving  laser  in  which  we  can  subsequently  drop  the  non- 
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secular  terms  with  respect  to  dressed  states.  In  a  previous  paper  we  executed 
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this  laborous  scheme,  and  it  turned  out  that  the  results  are  (almost)  identical 

to  the  simplified  derivation  which  we  now  give  here. 

In  order  to  achieve  an  expression  for  T  with  respect  to  dressed  states,  we 

merely  have  to  express  the  projectors  P  and  P  in  Eq.  (8.6)  in  dressed-states 

e  g 

basis  functions  |+>  and  |*>.  From  Eqs.  (10.5)  and  (10.6)  we  readily  find 


Pe  *  8-P+  ♦  8+P.  *  80(|+><-|  +  |-X+|)  (13.1) 

Pg  "  *+P+  +  +  8Q(|+X-|  +  |-X+|)  (13.2) 

where  we  have  introduced  the  abbreviations 
2  2 

g_  ■  cos  Jo  ,  8+  *  sin  J©  »  gQ  *  cosJOsinJ©  ,  (13.3) 

and  the  projectors  P  *  |t><±|  onto  the  dressed  states.  Next  we  substitute  the 

expansions  (13.1)  and  (13.2)  into  Eq.  (8.6),  and  drop  the  non-secular  terms. 
Care  should  be  exercised,  however,  because  a  J+>  -»  |+>  transition  between  two 
doublets  (Fig. 7)  is  in  exact  resonance  with  a  |->  -*•  |->  transition,  and 
therefore  couplings  between  these  transitions  should  be  retained.  Combining 
everything  then  results  in 

^  -  !  I 

kt 


+  g  va  -  a  +  a  -  a  )(P,pP  +  P  pP.) 
“o  eg  ee  ge  gg  +  -  -  + 


(13.4) 


where  the  term  proportional  to  gQ  comes  from  the  mentioned  degeneracy,  and 
denotes  a  projector  onto  a  dressed  state.  Representation  (13.4)  gives  rise  to  a 
master  equation  with  respect  to  dressed  states,  which  implies  that  we  can 
interpret  the  parameters  a^  as  the  rate  constants  for  transitions  between 


dressed  states. 
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Of  course,  the  parameters  a^  can  be  related  to  the  relaxation  constants 


with  respect  to  the  adbond  states.  For  k  *  e,g  and  t  *  e,g,  we  obviously 
have  a^  *  a^.  If  one  of  the  states  |k>,  |  i>  equals  a  |+>  or  |->  state,  we 
find 


*k,  ■  *j*k.  +  8t*kg 


(13.5) 


*tk  ’  »?*ek  +  *l*gk 


(13.6) 


in  terms  of  the  optical  parameters  g^  Transitions  between  the  |+>  and  |-> 
states  are  governed  by  the  rate  constants 


2  x  2  .  2/  ,  . 

g  a  +  g  a  +g(a  +a  ) 

t  ge  7  eg  °o  ee  gg 


it 


2,  ,  s  .  2  .  2 

g  (a  +  a  )  +  g  a  +  g  a 

Ro  ge  eg  8-  ee  gg 


(13.7) 

(13.8) 


Interaction  with  the  phonon  field  is  now  regarded  as  the  occurrence  of  single¬ 
phonon  transitions  between  dressed  states.  Relaxation  constants  which  connect 
an  adbond  state  jk>  *  |e>  or  |g>  with  a  dressed  state  |+>  or  |->  acquire  a 
contribution  from  two  distinct  processes,  due  to  the  fact  that  |k>  couples  with 
both  the  doublets  in  Fig.  7.  In  order  to  find  out  which  term  corresponds  to 
which  transition,  we  recall  the  definition  (7.5)  of  the  parameter  a^.  We 


notice  that  the  reservoir  correlation  function  5  is  evaluated  at  A,.,  which 

kl 


equals  the  frequency  of  the  phonon  for  that  particular  transition.  Therefore, 
the  terms  with  a^  and  a^e  represent  transitions  from  and  to  the  upper  doublet, 
respectively,  whereas  the  terms  with  ®gjt»aitg  describe  single-phonon  transitions 
between  |k>  and  the  lower  states.  In  a  similar  way  we  can  interpret  the  various 
terms  in  Eqs.  (13.7)  and  (13.8).  In  Fig.  8  the  different  processes  are 
indicated  by  arrows,  and  the  accompanying  optical  factor  specifies  the  term  in 
Eqs.  (13.5)  -  (13.8). 
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xrv.  PHOTON -PHONON  CONVERSION 

Stimulated  radiative  transitions  between  adbond  states  are  incorporated  in 
the  diagonalization  of  -  p-E,  and  therefore  only  single-phonon  transitions 
persist  in  a  pictorial  representation  of  the  various  processes  with  respect  to 
the  dressed  states.  In  order  to  elucidate  the  mechanism  of  photon  absorption, 
and  to  establish  the  relation  with  the  process  of  laser  heating,  we  consider  the 
licit  of  (relatively)  low  laser  power.  To  this  end,  we  first  note  that  g^  can 
be  expressed  in  ft  and  A  according  to 

g+  -  Hi  ?  (l  +  n2/*2)'*}  (I4.i) 

2  2 

which  depends  only  on  the  optical  parameters  through  the  combination  ft  / A  .  For 
weak  fields  we  then  obtain 

g_  -  1  -  fi2/4A2  ,  g+  -  ft2/4A2  (14.2) 

which  shows  that  g_  remains  present  without  a  radiation  field.  Rate  constants 

which  are  proportional  to  g_  must  consequently  correspond  to  radiationless 

transitions.  On  the  other  hand,  the  factor  g+  is  proportional  to  the  laser 

power,  which  implies  that  every  factor  g+  in  a  a^t  corresponds  to  the  absorption 

or  emission  of  a  photon.  In  this  fashion  we  can  track  down  the  significance  of 

the  optical  factors  in  the  relaxation  constants  with  respect  to  dressed  states. 

In  Fig.  9  we  draw  the  diagrams  with  respect  to  the  adbond  states,  and  they  are 

in  the  same  order  as  the  corresponding  diagrams  of  Fig.  8.  All  transitions  with 

2 

a  rate  constant  proportional  to  gQ*g  g+  are  single-photon  processes. 

Now  it  should  be  obvious  how  the  process  of  laser  heating  of  the  crystal 
can  be  conceived  as  a  result  of  many  photon-phonon  conversion  reactions.  Every 
phonon  in  the  diagram  corresponds  to  an  energy  exchange  between  the  adbond  and 
the  crystal,  whereas  a  photon  transition  amounts  to  an  energy  transfer  between 


the  laser  field  and  the  adbond.  A  combined  photon-phonon  diagram,  therefore, 
represent  effectively  an  exchange  of  energy  between  the  laser  and  the  crystal,  a 
process  which  is  mediated  by  the  adsorbate.  Summation  of  the  contributions  from 
all  diagrams,  weighted  with  the  probabilities  for  the  diagrams  to  occur,  then 
gives  the  net  power  absorption  by  the  crystal. 

XV.  SUMMARY 

We  have  studied  theoretically  the  optical  properties  of  an  adsorbed  atom, 
in  a  vibrational  bond,  on  a  crystal.  Coupling  of  the  adbond  states  with  the 
phonon  field  of  the  substrate  have  been  assumed  to  be  brought  about  by  single¬ 
phonon  transitions  only,  which  is  the  main  approximation  in  the  presented 
theory.  It  is  straightforward,  however,  to  include  higher-order  processes, 
especially  when  the  interaction  is  taken  to  be  a  Morse  potential.  Then  all 
matrix  elements  can  be  evaluated  analytically,  and  a  system  with  an  arbitrary 
number  of  adbond  states  can  be  parametrized  with  the  depth,  width  and  position 
of  the  binding  potential.  We  have  chosen  to  restrict  the  approach  to  single¬ 
phonon  transitions,  which  allows  a  clear  interpretation  of  the  coupling 
mechanism.  Then  we  can  incorporate  the  interaction  with  the  crystal  with  a 
single  phonon-field  amplitude  correlation  function  G(t),  which  was  studied  in 
detail  in  Section  IV.  It  is  the  behavior  of  this  crystal-response  function  that 
determines  whether  the  reservoir  approach  to  thermal  relaxation  can  be  justified 
or  not.  In  particular,  the  decay  of  G(x)  must  be  sufficiently  fast,  in 
comparison  with  the  adbond  relaxation  rates,  in  order  to  impose  a  Markovian  time 
evolution  on  the  reduced  adbond  density-operator.  This  decay  time  is  typically 
of  the  order  of  the  inverse  Debye  frequency,  which  is  reasonably  small  in 
comparison  with  the  rate  constants  for  a  single-phonon  process  (-  one  order  of 
magnitude) . 


Subsequently,  the  absorption  line  shape  has  been  evaluated,  showing 
considerable  deviation  from  a  Lorentzian  in  certain  cases.  The  origin  of  the 
distortions  and  the  shifts,  other  than  the  common  Lamb  shift,  has  been  tracked 
down  to  the  presence  of  non-secular  terms  in  the  time  evolution  of  the  density 
operator,  which  in  turn  appears  as  a  consequence  of  the  fact  that  the  damping 
constants  are  not  necessarily  small  in  comparison  with  the  transition 
frequencies. 

Next,  we  have  considered  the  irradiation  of  the  adsorbate  by  an  intense 
non -monochromatic  laser.  The  laser  linewidth  is  assumed  to  arise  from  a 
stochastically-fluctuating  diffusive  phase  of  the  driving  field.  Diagonialization 
of  the  Hamiltonian  is  accomplished  by  a  stochastic  transformation,  which  yields 
the  dressed  states.  These  states  can  be  interpreted  as  the  joint  eigenstates  of 
the  adbond,  the  single-mode  laser  and  their  interaction.  We  have  analyzed  the 
equation  of  motion  for  this  sytem  and  discussed  some  properties  of  the  steady 
state.  In  particular,  the  coherence  of  the  driven  transition  does  not  disappear 
in  the  long-time  limit,  which  indicates  that  the  system  is  not  in  thermal 
equilibrium. 

Then  the  power  absorption  has  been  obtained  from  the  work  done  on  the 
dipole  by  the  external  field.  The  finite  bandwidth  of  the  laser  gives  rise  to 


an  effective  photon  energy  with  ui^  in  between  and  Then  the  power 

absorption  can  be  written  as  I(w^)  ■  ^w^a^dl^-ll^) ,  with  a^  the  rate  of  optical 

transitions.  It  has  been  shown  that  is  linear  in  the  laser  intensity,  which 

implies  that  II  — H  tends  to  zero  in  the  high- intensity  limit,  since  the 
g  e 

absorption  rate  must  remain  finite.  Stimulated  transl  ions,  which  occur  at  a 


rate  a^,  Eq.  (11.4),  are  reduced  by  an  increasing  linewidth,  if  the  system  is 
driven  close  to  resonance  (6-0),  whereas  far  off  resonance  the  linewidth 
enhances  the  number  of  transitions  (a  -A).  This  feature  is  easily  understood: 


3 


in  the  case  A  -  0  most  laser  photons  are  farther  off  resonance  than  for 
monochromatic  incident  radiation.  On  the  other  hand,  for  a  large  detuning  there 
is  still  a  considerable  amount  of  photons  in  close  resonance,  which  have  a  large 
probability  to  be  absorbed. 

As  a  last  issue,  we  have  derived  the  rate  constants  for  single-phonon 
transitions  between  dressed  states.  In  the  low- intensity  limit,  these  processes 
can  be  interpreted  as  single-phonon/multiphoton  processes  in  adbond-state 
diagrams,  where  "multi"  stands  for  zero,  one  or  two.  This  reveals  that 
photoabsorption  is  inevitably  accompanied  by  a  downward  phonon  transition,  which 
gives  rise  to  heating  of  the  crystal.  Light- induced  desorption  through  resonant 
excitation  of  a  high-lying  state  is  consequently  bound  to  have  an  extremely  low 
ef  f iciency. 
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FIGURE  CAPTIONS 

Fig .  1 .  Real  (a)  and  imaginary  (b)  parts  of  the  correlation  function  G(t)  , 
divided  by  and  as  a  function  of  (d^t.  From  Eq.  (4.11)  it  follows  that  the 

integral  over  curve  (a)  equals  zero. 

Fig.  2.  Stimulated  part  of  the  reservoir  correlation  function,  divided  by 
_2 

2Cuijj  kT /Vi,  and  as  a  function  of  iDpt.  Curve  (a)  corresponds  to  lfaip  *  O.OOlkT 
(high-temperature  limit),  and  for  curve  (b)  we  have  taken  tfuip  *  4kT  (low- 
temperature  region). 

Fig.  3.  Real  (a)  and  imaginary  (b)  parts  of  the  Fourier-Laplace  transform  of 

-2 

the  spontaneous  part  of  the  correlation  function,  divided  by  Cu>p  *  and  as  a 
function  of  w/uip.  The  real  part  vanishes  for  w  <  0  and  for  <D><Dp,  whereas  an 
imaginary  part  is  present  for  every  <d. 

3 

Fig.  4.  Same  as  Fig.  3,  but  now  for  the  stimulated  part.  Normalized  by  CkT/Jtap 
and  for  huip  *  3kT.  The  real  part  only  vanishes  for  |uj|>(Dp.  Singularities  in 
the  imaginary  part  appear  at  u  *  i  which  is  an  artefact  of  the  sharp  cut-off 

of  the  dispersion  relation  at  oj^.  Any  smooth  decay  for  id  >  uip,  but  still 
arbitrarily  steep,  will  result  in  a  finite  value  of  the  correlation  function  at 

(D  *  ±  (Dp. 

Fig.  5.  Absorption  line  as  a  function  of  (d^,  and  divided  by  XpB(n^  - 
Frequencies  are  in  units  of  (dq.  The  parameters  for  these  curves  are  Imn  *  0  and 
Ren  ■  0.5, 1,2  for  a,b,c,  respectively.  For  a  decreasing  Ren  (-  linewidth),  the 
profile  tends  to  a  Lorentzian  with  half-width  at  half-maximum  equal  to  Ren*  For 
Imn  *  0  the  maximum  is  always  situated  at  (u^  *  <dq. 

Fig-  6.  Same  as  Fig.  5  but  now  with  Ren  M  0.2  for  all  three  curves.  The 
imaginary  parts  of  n*  the  lamb  shift,  are  +0.3,  -0.5  and  -0.6  for  a,  b  and  c 
respectively.  For  Imn  small,  the  line  around  *  <u  only  shifts  over  a 
distance  Imn,  but  for  larger  values  the  line  shape  changes  dramatically.  The 


top  of  the  curve  is  always  at  wL  =  wq( J l+2Imn/wo( )S ,  which  is  independent  of 
Re''.  For  Imn  small  in  comparison  with  the  top  appears  to  +  Imp. 

Fig.  7.  Diagonalization  of  the  dressed-adbond  Hamiltonian.  The  diagram  on  the 
left-hand  side  represent  product  states  of  the  adbond  |e>, |g>, |k>,. . .  and  the 
free  evolving  laser  |n>,  with  n  *  0,1,2...  the  number  of  photons  in  the  mode. 
States  |n>  and  |n+l>  are  separated  by  the  photon  frequency  In  this  diagram 

we  have  taken  A  *  >  0.  Since  the  laser  is  assumed  to  be  almost  on 

resonance,  we  have  a  ladder  of  doublets  separated  by  the  detuning  A.  Then  the 
dipole  interaction  couples  the  states  which  form  a  doublet,  but  it  does  not 
couple  two  sets  of  states.  A  diagonalization  in  turn  gives  rise  to  the  diagram 
on  the  right-hand  side,  where  the  separation  in  the  new  doublets  equals  n' , 
which  is  always  larger  than  A. 

Fig .  8 .  Single-phonon  excitations  of  the  laser-adbond  system.  The  diagrams 

with  the  arrows  in  the  opposite  direction  (not  shown  here)  correspond  to  phonon 

emission  into  the  crystal.  The  three  diagrams  with  the  double  arrows  persist  in 

the  absence  of  the  laser  and  correspond  to  radiationless  transitions  between 

adbond  states.  Solid  arrows  have  rate  constants  proportional  to  g+,  so  that 

2 

they  represent  single-photon  processes.  For  strong  fields  the  g+  transition 

appears,  which  couples  the  two  doublets  by  a  two-photon  process.  The  right-most 

diagrams  give  transitions  in  a  single  doublet  via  a  single-photon  process. 

Their  rate  costants  are  proportional  to  the  diagonal  matrix  elements  of  the 

derivative  of  the  binding  potential,  which  are  small.  For  zero  temperature 

these  transitions  vanish  identically,  since  then  we  have  a  *  a  *  0. 

ee  gg 

Fig.  9.  Low- intensity  interpretation  of  the  diagrams  from  Fig.  8.  With  respect 
to  the  adbond  states,  a  single-phonon  process  goes  together  with  photon 
absorptions  and  emissions  in  such  a  way  that  the  resulting  diagram  is  energy 
conserving.  The  sequence  of  processes  in  a  single  diagram  can  only  give  rise  to 


but  intermediate  states  can 
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transitions  between  the  real  states  |e>,  |g>, |k>, . . . , 


be  virtual  in  this  picture.  The  diagrams  are  in  the  same  order  as  in  Fig.  8, 


and  for  the  phonon -emission  processes  we  simply  reverse  the  directions  of  all 


arrows. 
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